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OPTIMIZATION  OF  THE  ACOUSTIC  PERFORMANCE  OF  A 
THICK  PLATE  WITH  EMBEDDED  ACTUATORS 

1.  INTRODUCTION 

The  use  of  piezoelectric  and  other  types  of  actuators  embedded  in  structures  to  either 
control  the  structure,  improve  its  performance  or  to  reduce  acoustic  noise  continues  to 
grow.  The  modeling  (see  [1]  for  instance)  and  optimization  of  the  global  performance  of 
these  “smart  structures”  [2]  is  an  important  concern. 

In  this  work,  the  focus  will  be  on  a  simple  plate  with  embedded  piston  type  actuators 
subjected  to  an  incident  acoustic  wave.  The  details  of  the  actuators  themselves  are  not 
modeled.  They  are  assumed  to  be  rectangular  sources  [3].  In  addition,  the  actuators  arc 
assumed  to  cover  a  partial  area  of  the  plate  with  an  impediance  equal  to  that  of  the  plate. 
This  assumption  could  be  relaxed.  The  present  work  relates  to  a  smart  panel  designed  for 
acoustic  quieting.  The  smart  panel  was  developed  as  part  of  the  Composite  Smart  Mater- 
rial^  (CSM)  Program,  which  is  an  interdisciplinary  effort  involving  several  private  corpo¬ 
rations,  academia  and  the  Naval  Research  Laboratory  [4]. 

Of  primary  concern  in  this  work  is  the  optimization/minimization  of  the  acoustic  in- 
tensityfor  pressure)  for  the  entire  far  field  as  well  as  in  selected  directions.  In  general, 
our  approach  is  limited  to  steady  state  conditions  because  of  the  availability  of  simple 
closed  form  solutions  for  scattered  and  radiated  pressures  at  specified  frequencies.  Tran¬ 
sient  solutions  for  this  class  of  problems  require  detailed  discretization  of  the  fluid  media 
which  is  computationally  very  costly,  and  lacks  accuracy  given  the  current  state-of-the- 
art.  In  fact  for  an  acoustic  fluid,  accurate  general  purpose  transient  boundary  integral 
numerical  solvers  are  presently  not  available  although  some  progress  has  been  made  re¬ 
cently  [5  ]. 

Optimization  is  itself  a  very  broad  topic.  Rao  [6]  presents  a  good  overview  of  optimi¬ 
zation  applied  to  enginnering  problems.  Unconstrained  optimization  techniques  have 
been  traditionally  divided  into  two  broad  catagories  -  direct  search  methods,  and  descent 
or  gradient  methods  [6,7].  Direct  search  methods  require  only  the  objective  function  val¬ 
ues  and  not  the  partial  derivatives  in  finding  the  minimum.  The  descent  methods  require 

derivatives  of  the  objective  function  and  are  known  as  gradient  methods. 
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Within  the  last  few  years  several  new  methods  have  emerged.  Among  these  are  neural 
network  techniques,  simulated  annealing  and  genetic  algorithms.  Neural  network  meth¬ 
ods  are  based  up  the  efficient  computing  power  of  the  network  of  interconnected  neuron 
processors.  Simulated  annealing  [6,7]  is  analogous  to  the  physical  process  of  annealing 
of  solids  and  has  been  successfully  applied  to  a  number  of  problems  among  them  the 
traveling  salesman  problem  [7].  It  is  well  suited  to  problems  of  large  scale,  especially 
where  a  global  minimum  is  hidden  among  many  poorer  solutions. 

Genetic  algorithms  (GA)  are  search  techniques  based  upon  natural  selection  and  ge¬ 
netics  [8-10].  They  can  be  quite  robost  and  have  been  successfully  applied  to  a  number 
of  difficult  problems.  We  will  be  applying  the  GA  to  the  problem  considered  here  and 
will  have  more  discussion  of  these  methods  in  section  2.2. 


2.  FORMULATION 

Consider  the  simple  plate  shown  in  Figure  1 .  The  plate  has  dimensions  2L,  by  2Ly, 
and  several  rectangular  piston  actuators  are  embedded.  The  dimensions  of  actuator  m  are 
by  bm-  Figure  1  indicates  4  of  these  macro-type  actuators,  which  are  comprised  of 
groups  of  smaller  pistons  [4],  but  any  reasonable  number  could  be  included.  In  addition, 
the  actuators  are  not  all  necessarily  at  100%  capacity,  nor  do  they  necessarily  provide 
coverage  of  the  entire  plate. 

As  indicated,  plane  wave  of  pressure  impinges  on  the  rectangular  plate.  Pjn^  is  of 
the  form; 


P.nc  = 


(1) 


where  i  -  Vm  ,  t  is  the  time,  O)  is  the  frequency  in  rad/s.  Pine  is  the  pressure  ampli¬ 
tude,  r  is  the  distance  and  k  is  the  wave  number  defined  by: 


k  =  0)/c  (2) 

In  equation  (2),  c  is  the  speed  of  sound  in  the  media. 

The  incident  wave  produces  a  scattered  wave  P^  at  a  field  point  which  in  spherical 
coordinates  ( r,6,^ )  is  expressed  as  [3]; 


-i2k  If  - 

Pjfr,©,^)  =  ^  e'  "^PjLjjLy  cos^j  •S(kLj^  •  (sin0j  cos^j  sin  0 sin 0)) 


(3) 
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(4) 


•S(kLy  (sinSj  cos^j  +sin0sin0)) 
where  the  function  S  has  the  form: 

0 

and  and  0j  represent  the  orientation  in  spherical  coordinates  of  the  incident  plane  wave. 

Referring  to  Figure  1  again,  consider  a  typical  actuator  m  with  dimensions  ajjjby 
b^and  its  centroid  locationed  at  (  ,  yom )  within  the  plate.  A  local  spherical  coordi¬ 

nate  system  (rm.^m.^m)  located  at  (xom.yom)-  From  reference  [3  ],  the  radiated 
pressure  Pfad.m  produced  by  actuator  m  being  driven  at  a  complex  velocity 

Vm=Vme"'“*  (5) 

is: 

Frad,m(rm,0in.^n,)  “  ‘  Prad,m(rm»0m>0ni) 

In  equations  (5)  and  (6),  Vm  is  the  velocity  amplitude  and  Prad.m  is  defined  as: 

Prad,m  (rm ,  0  m ,  )  =  -2iQ)p  V  m  •  (a^j  /  2)  •  (b,n  /  2)  •  S(.5ka,n  sin  0  m  cos^jj, )  (7) 

S(.5kbn,  sin  0 m  sin  ^  /(4;rr m  ) 

In  the  calculations,  a  field  point  is  expressed  in  the  global  spherical  coordinates  (r,0,0 ) 
and  then  transformed  to  the  m  actuator  local  spherical  coordinates  to  determine  Prad.m . 

2. 1  A  Standard  Mimimization  Approach 

With  the  scattered  solution  P^  from  equation  (3)  and  the  radiated  solution  P^aj  mfrorn 
actuator  m  given  by  equations  (6-7),  the  total  pressure  at  a  typical  field  point  (r,0,0 )  can 
be  expressed  by: 

M 

Ptot(r.0.0)  =  Ps(r,0,^)+  X*"rad.m(r.0.^)  (8) 

m=l 

In  this  equation,  the  M  actuators  may  or  may  not  cover  the  entire  plate.  Assuming  for  the 

moment  that  all  the  actuators  are  being  driven  at  the  same  velocity  amplitude  v ,  equation 
(8)  can  be  rewritten  as: 
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(9) 


_  M 

~  Prad.in(r,0,0) 

m=l 

The  factor  varies  between  0  and  1  and  allows  for  the  case  of  actuator  m  operating  at 
less  than  100%  capacity. 

Assuming  there  are  Nf  field  points  of  intere.st  in  the  far  field,  equation  (9)  can  be 
written  in  vector  form  as; 


—tot  ?s  '  ?rad 


where 


M 


—rad  ~  rad, 

m=l 


m 


(10) 


(11) 


In  equation  (1 1),  is  a  matrix  that  allows  the  field  point  locations  to  be  weighted  if  de¬ 
sired. 

We  next  define  a  performance  scalar  Jpthat  is  proportional  to  the  summation  of  the 
acoustic  intensities  of  the  far  field  [3]  as: 

.T 


P  ~  £tot  *  -tot 


where  implies  the  transpo.se  of  the  complex  conjugate  of  the  vector 
ducing  equations  (10-11)  into  (12)  produces: 

T  a  "P 

A*  _  —  _9=  = 

Jp=Ps  -Es+V  Prad*P,+V.pJ.P,3d-hV  P,ad-Prad 

From  [11],  the  minimum  value  of  v  is  determined  by: 

aT  aT 


(12) 

Intro- 


(13) 


'"’ir'  -"^rad* 


(14) 


•P,s)/(Prad*Prad) 

The  real  part  of  equation  (14)  reprc.scnts  the  velocity  amplitude  for  all  the  actuators  (some 
may  not  be  at  100%)  to  minimize  the  summation  of  the  acoustic  intensities  at  the  speci¬ 
fied  number  of  far  field  parts.  In  this  work,  we  will  refer  to  Vmin  produced  by  equation 
(14)  as  the  closed  form  minimum  (CFM)  velocity. 
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It  has  been  our  experience  that  if  there  is  a  significant  imaginary  component  of  Vmin , 
the  solution  is  quite  poor.  This  indicates  that  an  optimized  solution  with  all  the  actuators 
having  the  same  velocity  amplitudes  may  not  possible,  or  at  the  very  least,  require  the 
introduction  of  some  phase  differences  in  the  actuators. 

The  equations  in  this  section  were  programed  and  solved  using  the  commercially 
available  MATLAB  program  [12-13].  In  addition  to  selected  field  points,  our  resulting 
code  can  also  include  the  acoustic  intensity  of  the  entire  far  field  at  a  specified  radius  r. 
The  program  sets  up  in  MATLAB  a  grid  of  Nf  j  field  points  at  r  which  covers  the  limits 
of  6  and  <f) .  The  Jp  calculation  can  be  then  put  in  the  form: 

Jp  =  A.  Jpi /Nf +B  Jp2/Nfi  (15) 

where  A  and  B  vary  from  0  to  1  and  sum  to  1;  and  Jpi  and  Jp2  represent  the  summations 

that  are  proportional  to  the  intensities  of  all  specified  field  points  Nf  and  from  the  grid  of 
Nfj  field  points  covering  the  entire  far  field,  respectively. 

Returning  to  equations  (9)  zuid  (11),  one  can  allow  all  the  individual  actuators  to  be 
driven  at  different  amplitudes  Vm .  Equation  (9)  then  becomes: 

M  _  _ 

Pm  =  P,  +  Sf  in^mPrad,m  (16) 

m=l 

For  Nf  field  points,  equation  (16)  can  be  written  in  a  form  similar  to  (10)  or 


-tot  -s+^Srad,M 


(17) 


where  v  is  now  the  vector  of  the  actuator  velocity  amplitudes,  and  m  ^  matrix  of 

the  radiated  pressures  with  dimensions  Nf  by  M.  With  (17) ,  an  equation  similar  to  (12) 
can  be  generated  for  the  performance  scalar  Jp .  Application  of  standard  gradient  type 

techniques  [6,7,14]  can  then  be  introduced  to  solve  for  the  optimized  minimum  of  the 
actuator  velocities.  Our  experience  with  standard  gradient  techniques,  however,  proved 
very  disappointing.  The  solutions  became  stuck  in  local  minimums  unable  to  explore 
possibly  better  solutions  over  a  much  larger  solution  space. 
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2.2  A  Genetic  Algorithm 

Genetic  algorithms  (GA’s)  represent  a  type  of  stochastic  search  approach  that  attempts 
to  mimic  nature.  GA’s  are  well  suited  for  the  solution  of  many  optimum  design  problems 
that  are  characterized  by  mixed  continuou.s/discrete  variables,  and  large  discontinuous 
and  nonconvex  solution  spaces  [6].  Expensive  gradient  calculations  are  not  needed  for 
GA’s,  but  the  evaluation  of  individual  performances  must  be  inexpensive  and  fast  for  a 
genetic  algorithm  to  be  effective.  For  the  class  of  problem  we  are  investigating  in  this 
work,  these  above  conditions  are  satisfied. 

GA’s  are  based  upon  the  principles  of  genetics  and  natural  selection.  Basic  ingredi¬ 
ents  in  these  types  of  numerical  algorithms  include  reproduction,  crossover  and  mutation 
capabilities. 

For  the  problem  considered  in  this  work,  the  design  variables  will  be  the  velocity  am¬ 
plitudes  Vmof  the  indivdual  actuators.  The  velocity  amplitudes  are  transformed  into 
strings  of  binary  numbers,  0  or  1.  We  assume  here  that  each  string  is  6  bits  long  allowing 
2^  =  64  possible  discrete  values  for  each  actuator.  In  our  MATLAB  generated  program, 
both  lower  and  upper  velocity  amplitudes,  veil  and  vel2,  are  specified  for  all  the  actua¬ 
tors.  The  range  of  velocities  between  and  including  veil  and  vel2  is  divided  into  64  slots 
and  encoded  as  binary  strings.  Figure  2  indicates  an  example  having  4  actuators.  A  6  bit 
binary  string  for  each  actuator  velocity  is  shown  along  with  its  translation  into  the  appro¬ 
priate  velocity  slot.  This  24  bit  string  of  four,  6  bit  strings  represents  an  individual  in  the 
total  population. 

Reproduction  is  a  process  in  which  individuals  are  selected  to  be  parents  ba.se  upon 
their  fitness  measure  (such  as  Jp  in  equation  (13))  relative  to  the  rest  of  the  population. 

Each  individual  24  bit  string  j  is  assigned  a  probability  of  being  selected  for  reproduction 
via 

N 

Probj=:Jpj/^Jpjj  (18) 

ji=i 

where  Jp  j  is  the  performance  or  fitness  of  each  individual  and  N  is  the  size  of  the  popu¬ 
lation.  Highly  fit  individuals  have  a  greater  chance  for  reproduction,  while  less  fit  ones 
tend  toward  extinction. 
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Parent  1,  xi,  and  parent  2,  X2.  will  produce  two  offspring  yi  and  yj.  A  crossover  site  is 
then  selected  at  random  along  the  string  lengths,  and  the  binary  digits  are  swapped  be¬ 
tween  the  two  strings  following  the  crossover  site.  The  new  strings  yi  and  yi  are  then 
placed  in  the  new  population  and  the  process  continued  until  N  individuals  are  generated 
for  the  next  generation. 

A  mutation  operator  is  also  included  for  the  new  string  in  which  there  is  an  occasional 
random  alteration  of  a  binary  digit  in  the  string  from  a  0  to  a  1,  or  1  to  0.  This  is  con¬ 
trolled  by  a  mutation  probability  that  is  usually  set  low  (a  few  percent).  Mutation  is  a 
needed  ingredient  because  even  though  reproduction  and  crossover  effectively  search  and 
recombine  strings,  occasionally  some  important  genetic  information  can  be  lost,  or  the 
solution  can  be  trapped  at  local  minimums.  The  mutation  operator  helps  to  prevent  this, 
and  also  introduces  new  possibilities,  some  of  which  can  be  improvements,  into  the 
population. 

For  the  GA,  a  performance  scalar  Jp  defined  by  equation  (12)  rather  than  (15)  was  em¬ 
ployed.  Nf  selected  field  points  were  again  specified  for  the  far  field,  but  a  Jp2  term  was 
not  included  due  to  the  excessive  computational  cost  that  would  be  required  -  every 
member  of  the  population  in  every  generation  would  require  the  calculation  of  the  Jpi 
term. 

Another  point  to  be  mentioned  concerns  the  scaling  of  the  performance  (fitness) 
Jp  j  for  each  individual  in  the  population.  Outstanding  individuals  can  tend  to  dominate 

then  saturate  a  population  early  and  prevent  the  evolution  of  more  fit  individuals.  If  this 
occurs,  the  solutions  become  stuck  at  local  minimums  and  are  unable  to  more  fully  ex¬ 
plore  a  wider  range  of  the  solution  space.  One  simple  way  to  help  alleviate  this  situation 
is  to  modify  Jp  j  by  raising  it  to  a  power  Xn,  which  typically  is  in  the  range 

0+  <  Xn  <  1.0 .  The  quantity  Xp  varies  depending  upon  the  generation  number  n.  For  the 
early  generations,  x„is  set  to  some  specified  lower  limit  and  scaled  to  the  upper  limit  for 
the  last  generation.  Later  in  an  optimization  run  x^  approaches  the  upper  limit,  and  this 
allows  the  more  fit  individuals  to  then  saturate  the  population. 
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Finally,  we  note  that  because  genetic  algorithms  are  based  upon  probability  many 
times  it  is  necessary  to  make  several  computer  runs  to  obtain  best  results.  In  the  follow¬ 
ing  section,  this  is  implicit  in  the  discussion  of  the  results. 


3.  RESULTS  AND  DISCUSSION 

Consider  the  simple  plate  with  4  symmetrically  placed  actuators  as  shown  in  Figure  1. 
(The  numbering  of  the  actuators  begans  at  one  for  the  first  (Quadrant  then  proceeds  in  a 

counterclockwise  or  (p  direction).  An  incident  wave  Pjf|(,in  the  form  of  equation  (1) 

impinges  on  the  square  plate.  The  following  values  were  as.sumed;  Pjnc  =10  Pa, 
0inc  =7r/6  rad,  =0rad,  and  to  =  500  rad/s  for  the  incident  wave;  Lx  =Ly  =  1  m, 
a=b=.3  m  for  the  plate  and  actuators;  and  c  =  1500  m/s,  and  p  =  1000  kg/m’  for  the  fluid 
media.  Figures  3  and  4  are  intensity  Jp  (given  by  equation  (12))  plots  for  the  entire  far 
field  at  a  distance  of  r=500  m  for  the  closed  form  minimum  (CFM)  solution  from  equa¬ 
tion  (14);  and  the  GA.  For  both  solutions,  the  minimization  calculations  were  executed 
with  1^500  m  and  with  4  field  points  9=7tl6  and  <p  =  (0,7tl2,7tl6,nl2).  Figure  5  is  a 
baseline  plot  of  the  far  field  with  all  of  the  actuators  off. 

The  CFM  solution,  which  is  shown  in  Figure  3,  calculated  a  minimum  actuator  veloc¬ 
ity  of  5.17x10  ’  m/s  and  yielded  decibel  (dB)  reductions  of  the  intensities  of  the  entire  far 
field  and  4  field  points  of —7.4  dB,  and  —3.34,  -36.2,  -10.3,  and  —36.2  dB,  respectively. 
The  decibel  calculations  are  determined  from  the  formula  10*log(I/I,^f ),  where  I  and  I^r 
represent  the  acoustic  intensities  with  the  actuators  on  and  off  [15]  .  Many  other  combi¬ 
nations  of  various  field  point  locations  were  tried  as  well  as  including  the  entire  far  field 
via  equation  (15),  and  selecting  different  values  of  A  and  B.  The  results  in  Figure  4  arc 
representative  of  the  best  overall  minimization  of  the  far  field  by  the  CFM.  In  general, 
the  selection  of  far  field  points  in  direction  other  than  6  ^nl6  improved  tho.se  particular 
directions  only  so  much  with  minimum  impact  on  the  reduction  of  the  entire  far  field 
acoustic  intensity. 

Figure  4  indicates  about  the  best  results  obtained  via  the  genetic  algorithm  for  r=500 
m  with  4  field  points  6  =  n !  6  and  (f)  ~  (0,7r  I  2,k  /  6,71  /  2) .  The  population  size  for  the 
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GA  was  200,  the  number  of  generations  were  70,  and  a  fairly  high  mutation  factor  of 
0.10  was  used.  The  predicted  actuator  velocites  obtained  were:  .794x10'^,  9.05x10'^, 
8.41x10'^,  and  1.27X10'^  m/s.  The  dB  reductions  were  -3.79  for  the  entire  far  field,  and 
-4.6,  -7. 1 ,  -9.3  and  -7.2  dB  for  the  4  field  points.  In  obtaining  these  results,  the  GA  was 
limited  to  actuator  velocities  between  0  and  10.0x10'^  m/s.  Limiting  to  the  GA  to  a  nar¬ 
rower  range  closer  to  the  CFM  results  of  approximating  5.0x10'*  m/s  produced  GA  solu¬ 
tions  closer  and  closer  to  this  optimum  number. 

Overall,  based  on  these  and  numerous  other  results  for  both  the  CFM  and  the  GA,  we 
conclude  that  firing  all  the  4  actuators  at  approximately  5.0x10'*  m/s  produced  the  best 
overall  solution  for  this  symmetric  problem. 

We  next  consider  the  same  plate  and  incident  loading  with  again  4  actuators,  but  this 
time  the  actuators  were  running  at  100,  75,  10  and  50%  efficiencies,  respectively.  The 
far  field  intensities  Jp  were  optimized/minimized  for  8  field  points  or:  (0,0 )  =  (k/2,0), 

(7c/2,jc/2),  (7c/2,jt),  (jc/2, 1.5jc),  (7t/6,0),  (ji/6,jc/2),  (jc/6,7t/2)  and  (7c/6,1.57t).  Figures  6  and  7 
indicate  the  CFM  and  GA  Jp  (equation  (12))  results.  Although  quite  dissimilar  in  shape, 
the  CFM  and  GA  results  are  very  close  in  the  overall  reduction  of  the  intensity  of  the  en¬ 
tire  far  field  .  The  dB  reductions  were  -4.48  and  -4.45,  respectively.  For  these  GA  re¬ 
sults,  the  population  size  was  400  and  there  were  100  generations,  along  with  a  mutation 
factor  of  5%.  The  predicted  actuator  velocities  for  the  GA  were  7.62x10  * ,4.38x10'*, 
9.90x10'*,  and  1 1.81x10'*  m/s.  The  CFM  predicted  a  uniform  velocity  of  9.66x10  *  m/s 
with  a  significant  imaginary  component  of  .533x10'*  m/s. 

Actuators  2  and  3  (the  second  and  third  quadrants)  were  then  completely  turned  off, 
producing  an  even  more  unsymmetric  arrangement.  The  CFM  results  contained  an 
imaginary  component  in  the  velocity  of  actuators  1  and  4  that  was  approximately  25%  of 
the  real  part.  This  represented  a  very  poor  solution  that  actually  increased  the  intensity 
significantly  in  the  entire  far  field. 

Figure  8  shows  the  results  from  the  GA  with  actuators  2  and  3  turned  off.  The  GA 
produced  a  dB  reduction  of  the  entire  far  field  of  -3.4  dB  with  actuator  velocities  of 
7.16x10'*  and  7.00x10'*  m/s,  respectively. 

Using  only  the  genetic  algorithm,  we  next  consider  the  situation  of  16  symmetrically 
placed  actuators,  of  the  same  dimensions  as  before  with  a=b=  .3  m,  embedded  in  the 
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plate.  To  keep  the  actuator  velocities  similar  in  magnitude  to  the  previous  4  actuator 
models,  the  magnitude  of  the  incident  wave  Pjnc  was  increased  from  10  to  40  Pa,  while 
its  direction  remained  the  same  0,^^  =  rad  and  =  Orad,  and  to  =  500  rad/s.  The 
far  field  intensities  were  optimized  at  4  field  points  0  =  ;r /6  and  0  =  (0,n:/2,;r/6,;r /2) . 

The  population  size  was  200  and  70  generations  were  run.  The  GA  produced  actuator 
velocities  of  (1.43,3.49,3.97,  5.87,  2.38,  8.10,  8.25,  8.73,  8.89,  7.94,  9.84,  7.62,  6.51, 
3.97,  1.43,  10.0)xl0'^  m/s.  Figure  9  is  a  plot  of  the  entire  far  field  at  r=500  m,  and  Fig¬ 
ure  10  is  a  plot  with  all  the  actuators  off.  The  decibel  reductions  were  -6.42  dB  for  the 
entire  far  field,  and  -t-.178  (increase),  -13.76,  -14.8,  and  -13.7  dB  for  the  selected  field 
points.  Comparing  these  results  to  those  of  the  4  actuator  model  of  Figure  4,  we  see  that 
the  shape  of  the  far  field  has  changed  considerably  and  the  dB  reductions  are  greater  than 
previous  at  every  point  other  than  the  first  field  point  (0,0  )=(;r /6,0). 

Finally,  we  investigate  the  case  of  16  actuators  on  the  plate,  but  with  numbers  9-16 
turned  off  completely  (the  entire  second  and  third  quadrants).  In  an  attempt  to  cover 
more  of  the  far  field,  24  field  points  were  selected  with  0  and  0  varying  from 

e  =^inl2,nlA,7tl6)  and  0  =  (0,;r /4,;/r/2,.75;r,;r,1.25;r,1.5;r,1.75;r) .  Following  sev¬ 
eral  runs,  the  best  optimized  velocities  produced  by  the  GA  for  the  actuators  1  -  8  were 
(1.67,  .71,  1.67,  6.67,  5.95,  13.81,  10.95,  13.57)  xlO^  m/s.  Figure  1 1  indicates  the  entire 
far  field  acoustic  intensity  at  r=500  m.  The  decibel  reduction  for  the  entire  far  field  was 
computed  as  -4.7  dB,  which  was  actually  quite  good  considering  half  the  actuators  are 
turned  off  Note,  however,  that  as  in  Figure  9,  along  0  =  0  the  intensities  Jp  have  in¬ 
creased  slightly. 


4.  CONCLUSIONS 

In  this  work,  we  have  considered  a  simple  thick  plate  with  embedded  piston  type  ac¬ 
tuators  that  is  impinged  upon  by  an  acoustic  incident  wave.  The  actuators  do  not  neces¬ 
sarily  cover  the  entire  plate  nor  operate  at  100%.  The  question  addres.scd  concerns  the 
optimization  /minimization  of  the  acoustic  intensity  of  the  far  field  -  either  in  specific 
directions  or  the  entire  far  field. 
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A  standard  minimization  approach  was  first  employed  to  address  this  problem.  A 
closed  form  minimum  (CFM)  was  obtained  for  the  situation  in  which  the  actuators  are  all 
being  driven  at  the  same  velocity  amplitude.  The  reductions  in  the  far  field  were  excel¬ 
lent  for  the  situation  where  most  of  the  plate  is  covered  and  the  actuators  were  all  running 
at  near  100%  capacity.  For  even  some  moderate  non-symmetry  in  the  problem,  such  as 
indicated  in  the  discussion  for  Figure  6,  the  CFM  produced  reasonable  results.  However, 
the  CFM  results  declined  substantially  producing  erroneous  solutions  when  significant 
non-symmetry  was  introduced  into  the  problems.  In  these  cases,  equation  (14)  generated 

results  for  Vmin  with  substantial  imaginary  components.  This  indicated  that  an  optimized 
solution  with  all  the  actuator  velocity  amplitudes  the  same  may  not  be  feasible,  or  at  the 
very  least,  requires  the  use  of  phase  differences  in  the  individual  actuator  responses. 

The  use  of  standard  gradient  type  of  techniques  for  the  case  of  the  actuators  be  driven 
at  different  velocity  amplitudes  proved  to  be  ineffective.  Solutions  became  trapped  at 
local  minimums  unable  to  explore  larger  solution  spaces. 

The  genetic  algorithm  was  successfully  applied  to  all  the  cases  considered  in  this 
work.  The  GA  produced  feasible  solutions  in  which  all  the  actuators  were  driven  at  dif¬ 
ferent  velocity  amplitudes  (although  groups  of  actuators  could  have  been  tied  together  if 
desired).  For  a  symmetric  situation  the  GA  produced  smaller  dB  reductions  than  the 
CFM.  However,  unlike  the  CFM,  the  GA  generated  viable  solutions  for  the  reduction  of 
the  far  field  acoustic  intensities  in  all  cases  tried,  including  the  most  unsymmetric  situa¬ 
tions. 

Results  for  both  the  CFM  and  GA  approaches  revealed  that  only  limited  improvement 
in  the  far  field  intensity  reduction  was  possible  for  either  specific  directions  or  the  entire 
far  field.  In  addition  to  the  results  reported  in  this  work  in  which  and  were  fixed 

at  It/ 6  and  0  radians,  respectively,  numerous  other  values  and  combinations  of  actuators 
were  examined.  Those  results  tended  to  agree  with  this  statement.  However,  by  no 
means  did  we  conduct  an  exhaustive  numerical  study. 

In  the  future,  we  plan  on  seeking  to  include  delays  in  the  actuator  responses  in  both 
the  CFM  and  GA  form  of  the  equations.  This  should  improve  the  results  especially  for 
highly  non-symmetric  situations,  and  allow  further  reductions  in  specific  directions  by 
steering  the  response  into  other  directions  [15]. 
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6)  Figure  6.  CFM  intensity  results  for  4  actuators  running  at  100,  75,  10  and  50%  effi¬ 
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7)  Figure  7.  GA  results  for  4  actuators  running  at  100,  75,  10  and  50%  efficiency. 
Decibel  reduction  =  -4.45  for  the  total  intensity  of  the  entire  far  field. 
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8)  Figure  8.  GA  results  for  4  actuators  with  numbers  2  and  3  (second  and  third  quad¬ 
rants)  completely  off.  Decibel  reduction  =  -3.4  for  the  entire  far  field. 


9)  Figure  9.  GA  results  for  16  symmetrically  placed  actuators.  Decibel  reduction  =  - 
6.42  for  the  entire  far  field. 

10)  Figure  10.  Ba.seline,  all  1 6  actuators  are  off. 

1 1)  Figure  11.  GA  results  for  16  actuators  with  numbers  9-16  (second  and  third  quad¬ 
rants)  completely  off.  Decibel  reduction  =  -4.7  for  the  entire  far  field. 


14 


Incident 


Figure  1.  Simple  plate  with  embedded  actuators. 
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•Actuator  with  variable  velocities  (acoustic  solution,  same  frequency) 
•  Each  actuator  up  to  2^  =  64  possible  velocities  -  6  bits 
Actuator  Number  (assume  4  actuators  in  problem) 

1  2  3  4 

010100  111001  001001  101010 

-  20  =55  =9  =42  Example  values 


Figure  2. 
the  GA. 


Example  of  the  conversion  of  individual  actuator  velocities  to  binary  strings  for 
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Figure  3.  CFM  intensity  results  for  4  symmetrically  placed  actuators.  Decibel  reduc¬ 
tions:  total  intensity  of  far  field  =  -7.4;  selected  field  points  =  -4.6,  -36.2,  -10.3,  and  -36.2 
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Theta  (rad)  0  0  Phi  (rad) 


Figure  4.  GA  intensity  results  for  4  symmetrieally  placed  actuators.  Decibel  reductions; 
total  intensity  of  far  field  =  -3.8;  selected  field  points  =  -4.6,  -7.1,  -9.3,  and  -7.2  . 
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Figure  5.  Baseline,  all  4  actuators  are  off. 
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Figure  6.  CFM  intensity  results  for  4  actuators  running  at  100,  75,  10  and  50% 
efficiencies.  Decibel  reductions:  total  intensity  of  the  entire  far  field  =  -3.8;  selected  field 
points  =  -4.6,  -7.1, -9.3,  and -7.2. 
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Figure  7.  GA  results  for  4  actuators  running  at  100,  75,  10  and  50%  efficiency.  Decibel 
reduction  =  -4.45  for  the  total  intensity  of  the  entire  far  field. 
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Figure  8.  GA  results  for  4  actuators  with  numbers  2  and  3  (second  and  third  quadrants) 
completely  off.  DB  reduction  =  -3.4  for  the  entire  far  field. 
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Figure  10.  Baseline,  all  16  actuators  are  off. 
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Figure  11.  GA  results  for  16  actuators  with  numbers  9-16  (second  and  third  quadrants) 
completely  off.  Decibel  reduction  =  -4.7  for  the  entire  far  field. 
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